GRAVITATIONAL LENS STATISTICS FOR GENERALIZED 

NFW PROFILES: PARAMETER DEGENERACY AND 
IMPLICATIONS FOR SELF-INTERACTING COLD DARK 

MATTER 

J. S. B. Wyithe 1 ' 2 , E. L. Turner 1 & D. N. Spergel 1 ' 3 
Received ; accepted 



1 Princeton University Observatory, Peyton Hall, Princeton, NJ 08544, USA 
2 School of Physics, The University of Melbourne, Parkville, Vic, 3052, Australia 
3 The Institute for Advanced Study, Princeton, NJ 08540 



- 2 - 



ABSTRACT 

Strong lensing is a powerful probe of the distribution of matter in the cores 
of clusters of galaxies. Recent studies suggest that the cold dark matter model 
predicts cores that are denser than those observed in galaxies, groups and 
clusters. One possible resolution of the discrepancy is that the dark matter has 
strong interactions (SIDM), which leads to lower central densities. A generalized 
form of the Navarro, Frenk and White profile (Zhao profile) may be used to 
describe these halos. In this paper we examine gravitational lensing statistics 
for this class of model. 

The optical depth to multiple imaging is a very sensitive function of the 
profile parameters in the range of interest for SIDM halos around clusters of 
galaxies. Less concentrated profiles, which result from larger self-interaction 
cross-sections, can produce many fewer lensed pairs. Furthermore, profiles that 
result in a small optical depth exhibit reduced typical splittings, but produce 
multiple images that are more highly magnified. However the resulting increased 
magnification bias does not alter our conclusions. 

We find that lensing statistics based on profile parameters obtained from 
fits out to the virial radius are dependent on the minimization scheme adopted, 
and may be seriously in error. However, profile fits weighted towards the core 
region have parameter degeneracies that are approximately equivalent to those 
for strong lensing cross-sections. 

Lensing statistics provide a powerful test for SIDM. More realistic and 
observationally oriented calculations remain to be done, however larger 
self-interaction cross-sections may well be ruled out by the very existence of 
strong lenses on galaxy cluster scales. The inclusion of centrally dominant 
cluster galaxies should boost the cross-section to multiple imaging. However 
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our preliminary calculations suggest that the additional multiple imaging rate is 
small with respect to the differences in multiple imaging rate for different halo 
profiles. In future statistical studies, it will be important to properly account 
for the scatter among halo profiles since the optical depth to multiple imaging 
is dominated by the most concentrated members of a cluster population. 

Subject headings: clusters: general - gravitational lenses - dark matter - halos 
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1. INTRODUCTION 

The cold dark matter plus cosmological constant model (ACDM) appears to be 
remarkably successful in describing cosmology on scales much larger than a Mpc (e.g. 
Bahcall, Ostriker, Perlmutter & Steinhardt 1999). However recent observations on smaller 
scales find ACDM wanting: ACDM predicts dwarf galaxy halos (e.g. Moore et al. 1999) 
that are denser and more concentrated than observed (Dalcanton & Bernstein 2000; 
Firmani et al. 2000b). The persistence of bars in galaxies like our own, suggests that the 
dark matter is not centrally concentrated (Sellwood 2000). In addition, the intensively 
studied CL0024+1654 is nearly spherical (Colley, Tyson & Turner 1996; Tyson, Kochanski 
& Dell'antonio 1998), at odds with ACDM predictions of triaxial systems. Spergel & 
Steinhardt (2000) have recently suggested that these and other discrepancies (see for 
example Dave, Spergel, Steinhardt & Wandelt 2000; Wandelt et al. 2000 and references 
therein) could be resolved if cold dark matter were weakly self interacting with a large 
scattering cross-section (SIDM). Several authors (Burkert 2000; Yoshida, Springel, White 
& Tormen 2000; Dave et al. 2000) have since undertaken N-body simulations of halos with 
self interactions in the regime suggested by Spergel & Steinhardt (2000). These works find 
that SIDM halos produce flatter, smoother and more spherical cores than their ACDM 
counterparts. 

We are interested in the statistics of multiple imaging due to a generalized NFW profile 
(Zhao 1996) which may be used to describe both CDM and SIDM halos. We concentrate 
our efforts on lensing by clusters of galaxies, since single galaxies contain a significant mass 
in baryons, distributed to produce a flat rotation curve. This makes SIDM in galaxies 
difficult to probe using gravitational lensing. Conversely, clusters of galaxies have a much 
larger mass to light ratio, and baryons which are less centrally concentrated (particularly 
in the absence of a central dominant galaxy). We assume spherically symmetric profiles 
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and concentrate on the statistics of the optical depth to multiple imaging. The total 
magnification, and the maximum image splitting are also considered. 

In Sees. |2] and |3] we discuss the properties of dark matter halos and introduce the 
formalism for computing the lensing statistics of the Zhao profile. In Sees. U |5] and |7| we 
discuss the statistics obtained, contrast the parameter degeneracies present in the lens 
statistics and in profile fitting, and present some implications for SIDM cluster halos. In 
Sec. H we discuss the important implications for lens statistics of the distribution of halo 
profile-parameters. We also discuss two potentially important caveats which turn out not 
to effect our conclusions: we discuss magnification bias (Sec. [(]), and give a preliminary 
account of the effect of massive central galaxies on the multiple imaging rate (Sec. |]). 
Throughout the paper we have assumed a standard Q a — 1, A = and h = 100k ^° sec 1 =0.7 
filled beam cosmology. 



2. DARK MATTER PROFILES 

Navarro, Frenk & White (1997) claimed the NFW profile (with an inner density p(r) 
varying with radius r as p(r) oc r~ l ) to be universal. More recently Moore et al. (1999) 

3 

and Ghingna et al. (2000) have claimed a different universal profile, having p{r) oc r~2 in 
the innermost regions. However Jing & Suto (2000) find a systematic correlation of inner 
profile slope with mass. Moreover, they find some scatter in the slope at a fixed mass. On 
the other hand, from their convergence study, Klypin, Kravtsov, Bullock & Primack (2000) 
find that while there is some scatter, in general halos are equally well fit in the region 
where convergence is achieved by either an NFW or a Moore profile, independent of the 
halo mass. They interpret the correlation found by Jing & Suto (2000) as a trend in halo 
concentration rather than in asymptotic slope. Dave, Spergel, Steinhardt & Wandelt (2000) 
present a comprehensive suite of simulations of both CDM and SIDM dwarf galaxy halos. 
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They find typical profiles for CDM halos in the inner most regions of p(r) ~ r~2, with some 
scatter. However they found SIDM halos to be less cuspy, having central profiles ranging 
from p(r) ~ r^ 1 to p(r) ~ r~5. The uncertainty as to the optimal halo profile as well as 
the variation introduced by SIDM motivate an exploration of lensing by profiles that differ 
from the analytically treatable NFW form. 

Following the introduction of the Navarro-Frenk- White profile (NFW) (Navarro, 
Frenk & White 1995,1996,1997), Bartelmann (1996) was the first to consider its lensing 
implications, showing that the profile is capable of producing radial arcs like those observed 
in several clusters (e.g. Mellier, Fort & Kneib 1993; Smail et al. 1995). There have 
been many subsequent efforts contrasting resultant lensing predictions with observation. 
Bartelmann et al. (1997) found the optical depth for large arcs for lensing by NFW profiles 
and found significant dependence on cosmology, though no models could reproduce the 
observed number. Wright & Brainerd (2000) and Asano (2000) have contrasted weak lensing 
mass measurements and 2-image system image positions flux ratios based on the NFW 
with those for a singular isothermal sphere (SIS). Maoz, Rix, Gal- Yam & Gould (1997) 
carried out a survey for wide separation lensed quasars, and compared their results with 
a thorough statistical study of the number of expected cluster lensed quasars. Modeling 
clusters using the NFW profile, they found that their null result for observed lensed quasars 
was consistent with the lensing statistics. Molikawa, Hattori, Kneib & Yamashita (1999) 
investigated giant arc statistics and found that the NFW profile can reproduce the number 
of observed lensed arcs. However their study finds some discrepancy between virial and 
X-ray temperature for some clusters. Broadhurst, Huang, Frye & Ellis (2000) have found 
that an NFW profile (singular core) can explain strong lensing in CL0024+1654, contrary to 
earlier claims of a flat core (Tyson et al. 1998). However Shapiro & Iliev (2000) find that the 
velocity dispersion implied by the NFW fit is a factor of two larger than that observed. The 
shallower cores resulting from SIDM simulations may resolve this discrepancy. On the other 



-7- 



hand, Smith et al. (2000) find that radial arcs imply a core in A383 which has a steeper 
logarithmic slope than an NFW profile. Williams, Navarro & Bartelmann (1999) find that 
arc properties for the most massive clusters are consistent with NFW profiles derived from 
N-body simulations, but that clusters with velocity dispersions of ~ 1000 km sec' 1 require 
the presence of a massive central galaxy to reconcile the models with observation. Fox, 
& Pen (2000) have considered lensing by galaxy groups modeled as NFW profiles in the 
Hubble Deep Field. They find that the null result for lenses is consistent with lensing 
probability. In addition, they compare the probability for multiple imaging due to an NFW 
profile with the probability for an SIS, and find that the probability is significantly reduced. 

Recent complimentary studies have also considered lensing statistics for Zhao profiles. 
Molikawa & Hattori (2000) find that the ratio of radial to tangential arcs is a sensitive 
function of the inner slope of the density profile. Li & Ostriker (2000) consider a 
semi-analytic treatment of lensing by the Zhao profile. They find that lensing probability is 
sensitive to the profile considered, as well as cosmology. In addition, they find a probability 
for large splittings which is inconsistent with the observed distribution. 



3. GRAVITATIONAL LENSING BY GENERALIZED NFW PROFILES 

The NFW halo density profile can be generalized to describe a profile p{r) as a function 
of radius r with an arbitrary power-law shaped central cusp p(r) ~ r _/3 and outer regions 
that fall off as r~ 3 (Zhao 1996), 

fe) + 

Navarro, Frenk and White (1996) defined the concentration parameter 

Cnfw = — (2) 
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where 

is the radius enclosing a mass M 20 o with an over-density of 200 above the critical density. 
The characteristic over-density S c in Eqn. [I] is 

^ _ 200C^r FV1/ 

3F(CVfwO 

where 

F(y) = [ V x 2 -P(l + xf- 3 dx. (5) 
Jo 

Therefore, given a cluster mass M 20 o, the Zhao profile has two free parameters (3 and 
Cnfw- Cnfw is quoted for a profile at red-shift zero. Our lensing calculations assume 
a non-evolving cluster-lens population. The concentration parameter Cnfw{z) and 
characteristic over-density S c (z) of a halo placed at red-shift z are therefore the solution of 

C NFW {z) ■ C NFW (z)- f{Cnfw) {1 + z)2{1 + noZ) (6) 



and 



These transformations leave r s and p c S c independent of z. From the study of the evolution 
of a large sample of ACDM halos, Bullock et al. (2000) find C NFW {f3 = 1) oc (1 + z)~\ 
consistent with Eqn. |6|. Note however, that the inclusion of a finite self-interaction 
cross-section should result in a different profile evolution. 

The Zhao profile is spherically symmetric, while real clusters, and N-body simulations 
tend to be elliptical or tri-axial. Blandford & Kochanek (1987) and Kochanek & Blandford 
(1987) find that the the introduction of ellipticities smaller than 0.2 into nearly singular 
profiles has little effect on the cross-section and image magnification. Furthermore, the 
images remain on opposite sides of the potential so that the variation in image separations 
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is similar to the axis ratio of the potential. However profiles containing a finite core and 
ellipticity show a significant increase in cross-section and qualitative differences in image 
positions. 

We make the thin screen approximation for gravitational lensing by extended bodies. 
The surface mass density at a radius £ may be written 



S(0 = 25 cPc r s x l - f3 / 2 sin9(sm9 + xf~ 3 d9, 
Jo 



where x = The mass M(£) enclosed within the cylinder of radius £ and the bend angle 
a(£) of a ray with impact parameter £ are 

Af(£) = 2vrr s 2 [ x'Y,(x')dx' . (9) 

and 



o 



a(£) = (10) 



If (3 = 1, Eqns. ||, H and reduce to the expressions computed by Bartelmann (1996). 



The lens equation describes the location £ of an image in the lens plane given a source 
position rj (both measured from the observer source line of sight). The angular diameter 
distances from the observer to a lens at red-shift Zd, from the observer to a source at z s and 
from the lens to the source are denoted by Da, D s and Dd s . With these, the lens equation 
may be written (e.g. Schneider, Ehlers & Falco 1992) 

= 9-(3-a(6)^, (11) 

where 9 = -g- and j3 = are the angular positions of the image and source respectively. 
The point source magnification ^ of an image at 9{ is 

"■^(r^W 1 ' (12) 
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The Zhao profile has a critical radius r] crit separating regions of the source plane that 
produce 1 or 3 images. The corresponding critical impact parameter ^ crit is the solution of 



(13) 



and the critical angle for multiple imaging is 



^u = ^-a(0^f. (14) 



For later comparison, the bend angle for a SIS is 



aSIS = 2ttGJU™ 
c 2 r 200 

In this case = so that 

nSIS _ 2kGM 200 D ds 

Peril - 2 n ■ l i0 J 

c r 2 oo L) d 

For the computation of lens statistics for an unevolving population of lenses with 
constant co- moving number density n d = n (l + z) 3 we find (Turner, Ostriker & Gott 1984) 
that for an Q = 1 universe the differential and total cross-section may be written 

dr(z d , z 8 ) = im —D 2 d (z d )l3 2 crit (z d , z s )y/l + z d dz d (17) 

and 

vet/. D ^Mnt( Zd, z s )y/l + z d dz d . (18) 

ti Jo 

Turner, Ostriker & Gott (1984) defined a dimensionless quantity F S is that measures the 
efficiency of cosmologically distributed SISs for producing multiple images. 



"sis 



4^ (19) 



We use Fsis to normalize r. More realistic calculations of t(z s ) for the NFW {(3 = 1) profile 
have previously been made within the contexts of wide separation quasars (Maoz, Rix, 
Gal- Yam & Gould 1997), giant arc statistics (Bartelmann et al. 1997; Molikawa, Hattori, 
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Kneib & Yamashita 1999) and distant supernovae (Porciani & Madau 2000). However in 
this work we are interested in relative values of r produced by different profiles, and so have 
neglected quantitatively important effects such as formation history. 

One image is present for all source positions. We denote this image to be at 9\. Note 
that 6*i > P cr i t for all [3. For (3 < (3 crit there are two additional images with positions 
9 2 < (3 cr it and 6*3 > ficrit- We label the corresponding point source image magnifications 
as fii, /i 2 and ^3 respectively. If (3 crit > then 9\ > 0, 61, 9 2 < 0, fj,%, fi 3 > and /i 2 < 0. 
In addition \fi 3 \ > |/i 2 |- For (3 < (3 crit , we define the average image separation of a 
multipli-imaged source at z s due to a deflector at z d to be 

{A9(z d , z s )) = -2- / f3{9 x {(3, z d , z s ) - 9 3 (J3, z d , z s ))d(3. (20) 

Per it •' 

Similarly the average magnification of a multipli-imaged source at z s due to a deflector at 
Zd is 

2 [Pent 

(fi(z d ,z s )) = -3— / /3(|^i(/3,2: d ,0 s )| + |/i 2 (^,£d,^)| + |^ 3 (^,^,^)|)c?/3- (21) 

Pcrit J ° 

It follows that the average image splittings (A9(z s )) and total magnifications (n(z s )) for a 
source at z s are 

(A9(z s )) = [ Zs (A9(z d ,z s )) dT(z _ d > Zs) dz d (22) 

T{z s ) Jo dz d 

and 

<M*)> = "7~T l Z \^z s )) dT{Zd ' Zs) dz d . (23) 

4. CROSS-SECTION, IMAGE SPLITTINGS AND MAGNIFICATION 

In this section we present the differential and total cross-sections (Eqns. |H] and |IBD , 
distributions of image separations, and distributions of magnifications (Eqns. and |23|) for 
Zhao profiles having different parameters (3 and Cnfw- 
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Fig. [T] shows p(r) (top) and £(£) (bottom) for three different Zhao profiles. These are 
labeled profile / (solid lines), II (dashed lines) and III (dotted lines), have parameters of 
(/3, Cnfw) = (1-5,20), (1.0,15), and (0.5,10) respectively, and are increasingly steep cusps 
in their central regions in addition to being increasingly centrally concentrated. The model 
clusters have an M 2 oo equal to that of a SIS (M^qq) with a velocity dispersion corresponding 
to a large cluster cr'gjg = 1000 km sec -1 : 

^=<^K»y '(§)'• « 24 » 

The solid light lines show psis{ r ) an d ^sis(0 °f the SIS for comparison. The horizontal 
dashed lines show the critical surface mass density 

c 2 D s 

EcrU = A^GD^ (25) 

for a source at z s = 3 and lenses at Zd = 0.5, 1.0 and 1.5. Comparison of S cr j t with 

indicates a significant variation of the lensing probability between profiles having the 
same M 20 o but different profile parameters. Profiles / and II reach the critical density 
over a much larger area than profile III, which is only just supercritical at Zd = 1-0, and 
sub-critical at z,i = 2.0. 

The variation of cross-section with the profile parameters is quantified by Fig. ^| which 
shows t(z s ) for Zhao profiles / — III. The cross-section is given in units of F$is (Eqn. [19]). 
For comparison, the cross-section for a SIS (tsis( z s)> light line, Turner, Ostriker & Gott 
1984) is also shown in Fig. |[ Examples of the effect of different cosmologies on tsis{z s ) 
(and on 7(~^ ^dz'd^ m can ^ e f° un d i n Turner et al. (1984) and Turner (1990). 

There is a striking variation in r over these three profiles. Profiles II and III produce rs 
which differ by an order of magnitude, while there are 4 orders of magnitude between r 
for profiles I and More concentrated profiles are much more efficient for producing 
multiple images due to their higher central densities. For larger source red-shifts the SIS has 
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a similar cross-section to profile //, while at smaller red-shifts the cross-section is similar to 
that of profile /, indicating that different profiles are efficient lenses at different red-shifts. 

In Fig. |^ we have plotted dT ^ a z ^ for profiles / (solid lines), II (dashed lines) 
and III (dotted lines) with z s = 1.5 (upper panel) and z s = 3.0 (lower panel). The 
light lines show the corresponding distributions for a SIS. Less concentrated profiles have 
distributions that are more sharply peaked, with modes at lower red-shift (c.f. Kochanek 
& Blandford 1987). The lower cross-section to multiple imaging is therefore due to a 
combination two effects. Firstly, the halo reaches the critical density to lensing only over a 
small region in the center due to the low concentration or central slope. In addition, there 
is a decreased redshift range in which the central surface density exceeds the critical value 
anywhere. Note that the assumption in Eqns. |17] and [18] of an unevolving population of 
SIDM halos is aided in the case of low-concentration profiles since these preferentially lens 
at low red-shift. 

To further investigate the variation of cross-section with (3 and Cnfw we have plotted 
contours of T ^(l \ with z s = 3, over a parameter space of (3 and Cnfw (Fig. [|). Two 
additional dark lines are shown. The upper dark line is a contour of unity and marks the 
parameters that produce an cross-section equal to that of a SIS. The lower dark line marks 
the boundary of the region of parameter space where no multiple imaging is produced (on 
this and subsequent plots this line is the approximated by the 10~ 5 contour). The contours 
describe sets of parameters (3 and Cnfw that are equally efficient for producing multiple 
images. Profiles /, II and III lie on a line in f3 — Cnfw space that is approximately 
orthogonal to the contours of t(z s ). However, figure f| demonstrates that dramatic 
variability in t(z s ) may be obtained by variation of either (3 or Cnfw- 

We have calculated \ , (A0(z a )) and (n(z s )) over a small grid of (3 and Cnfw (see 

Tab. [I]). We find that profiles which produce a small cross-section to multiple imaging also 
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produce smaller image splittings, although only by a factor of a few compared with more 
highly concentrated profiles. However, the image magnification is significantly larger for the 
less concentrated profiles. Distributions of image separations ( d ^ z y ) for the Zhao profile 
are presented in Fig. || for the values of j3 and Cnfw m Tab. |I[ For comparison, the light 
line shows the distribution of splitting angles for a SIS (^w^y)- There is a trend for the 
distributions with lower (A9(z s )) to also have a narrower range of A8(z s ). In contrast to 
the SIS distribution, the Zhao profile distributions are asymmetric about the mode, having 
a longer large splitting tail. Fig. || shows magnification distributions ^ for the values of (3 
and Cnfw i n Tab. [l|. The distribution for an SIS, 

dP s,s 8 , , 

— = 7' (26) 

is shown for comparison (light line). From Fig. |]we see that the high mean magnifications 
shown in Tab. [I] result from a higher minimum magnification of multiple images. Lenses 
with low concentration and small f3 result in a low t(z s ). Very close alignment is therefore 
required to produce multiple imaging, and the close alignment results in high magnification. 
More concentrated lenses do not require such high alignment for multiple imaging. This 
results in a larger cross-section and reduced mean magnification. 

The basic lensing properties of the Zhao profile may therefore be summarized as follows. 
Profiles that have a small concentration and flat central profile produce significantly less 
multiple imaging than a profile with either a high concentration or steep central profile (or 
both). However, these few multiple images are significantly magnified. In any lens survey, 
the increased (fx) will partially offset the decreased r(z s ) in determining the lens frequency. 
This magnification bias is discussed in Sec. |6|. For clusters of galaxies, the reduced splittings 
for low density profiles are ^ 10" and should not pose a significant selection bias for well 
defined surveys. 
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5. THE C NFW - ft DEGENERACY 

A property of the Zhao profile (Eqn. [1]) is the ability to describe a high central density 
either by a steep central cusp, or by a large concentration parameter (so that the turnover 
separating the core regions from the outer halo occurs closer to the center). This degeneracy 
is present between other profiles exhibiting a power-law core and a characteristic radius. 
Klypin, Kravtsov, Bullock & Primack (2000) have investigated the relative suitability of 
Moore et al. (1999) and NFW (ft = 1) profiles as descriptions of CDM halos. For galaxy 
size halos, they find differences which are small at radii above 0.01r 20 o- Larger differences 
are found for lower concentration cluster halos as the profiles central slope begins to 
approach its asymptotic value. Moreover, using a very general 4-parameter NFW profile, 
Klypin et al. (2000) obtained equally good fits to an N-body halo for several different sets 
of parameters. Our interest lies in the relationship between the Cnfw — ft degeneracy found 
in in the contexts of multiple imaging statistics and profile fits to simulated halos. 

Our approach is to select a collection of profiles p (r) having ft = 0.5 but different 
concentrations Cnfw,o- Keeping M 2 oo constant (for comparison with our strong lensing 
calculations), we fit a new profile Pfit(f) with Cnfw as a function of ft. Fits are obtained 
by minimizing each of three different quantities which emphasize different parts of the halo: 

* 2= Lo r 2 (p(r) - Po (r)) 2 dr, (27) 

Too" 

kpmax = max(\log(p(r)) - log(p (r))|) (28) 
(Klypin, Kravtsov, Bullock & Primack 2000), and 

Minimization is performed over the range 0.01r2oo — r 2oo- 

Some examples of the fitted profiles and their fractional residuals are presented in 
Fig. |7| (light lines). The dark lines show the profile being fitted. In each case we quote the 
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parameters O , Cnfw,o and /3, Cnfw, as well as the average absolute relative difference over 
the fitting volume 



The central density dominates minimization of x 2 (top panels) so that the fitted profile 
has the correct central density at the cost of error over large regions around r = 
Alternatively if we set Cnfw, an d adjust /3, the fitted profile has the correct central slope 
and density, but is systematically high or low over the outer regions. The initial and fitted 
halos have average absolute relative differences of a few to ~ 15%. As noted by Klypin, 
Kravtsov, Bullock & Primack (2000), minimization of Ap max (central panels) often does a 
much better job of reproducing the density at small r. If (3 is the fitted parameter rather 
than Cnfwi a similar result is obtained, the central slope adjusts so as to give the correct 
density at small r. Minimization of Ap max produces halo fits with similar average absolute 
relative differences to x 2 . Finally, minimizing x 2 e i (lower panels) obtains good fits beyond 
r ~ 0.1r2oo at the expense of an incorrect central density. The mean relative error is 
substantially reduced, with average absolute relative differences of a fraction of a percent. 

The dot-dashed lines in the upper, central and lower plots of Fig. |] are depictions of 
the [3 — Cnfw degeneracies corresponding to the minimization quantities x 2 ■, Ap max and 
Xrei- Since x 2 produces good fits to the core region, which is known to be most important 
for strong lensing (e.g. Turner, Ostriker & Gott 1984), the lines of degeneracy are nearly 
parallel to those describing the cross-section to multiple imaging over much of the parameter 
space. However for combinations of low or high concentration and low inner slope, the 
approximately degenerate profiles do not provide a accurate description of the lensing 
statistics. Minimization of Ap max results in profile degeneracies similar to those for strong 
lensing unless the profile has large C^fw an d small [3. Finally, minimization of x 2 e i leads 
to degeneracies that differ significantly from those for strong lensing. Fits that are weighted 




(30) 
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towards the central regions, and which result in highly concentrated profiles may therefore 
be converted to the analytically treatable NFW {(3=1) profile while approximately 
preserving the lensing statistics. This is particularly true for fits based on minimization of 
^Pmax- On the other hand, lensing calculations based on profiles which best describe the 
density over the largest volume (out to r 2 oo) may severely under or over-estimate the strong 
lensing statistics. For example, Zhao halos having the same M 2 oo, but different parameters 
(P,Cnfw) — (1.5,6.0) and (1.0,10.0) are degenerate to fitting by the minimization of xleii 
but differ in r by an order of magnitude. Thus Fig. [8] may be considered a quantification 
(for Zhao profiles) of the statement that lensing is sensitive to the size of the halo core. 



6. MAGNIFICATION BIAS 

Statistics generated from surveys for gravitational lenses are subject to several classes 
of bias. Kockanek (1991) has given a detailed discussion of selection biases in optical, 
imaging, galaxy mass gravitational lens surveys. Two primary sources of bias against 
finding a lens are selection against small separations due to the finite angular resolution of 
the survey, and the obscuration or differential reddening of one or both images due to dust 
in the lens galaxy. When considering lensing by clusters, the splittings are nearly always 
larger than the survey resolution. In addition, the lensed images are viewed through the 
halo rather than the outskirts of a galaxy. These biases should not therefore operate on any 
significant level in the analysis of cluster lensing statistics. The magnification of a lensed 
image has the potential to boost the observed flux of a source which is intrinsically faint 
above the survey detection limit. Lensed images therefore sample fainter members of the 
source population. This is magnification bias (Gott & Gunn 1974; Turner 1980), and it a 
significant factor in any lens survey. 

As mentioned in Sec. f|, lenses with low concentration and a shallow cusp result in a 
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low cross-section. Very close alignment is therefore required to produce multiple imaging, 
and this results in a high minimum magnification. More concentrated lenses do not require 
such high alignment to produce multiple imaging, thus the cross section is larger, and 
the minimum magnification not so high. In the latter case, although the mean is lower, 
high magnifications are still present where the alignment is high. As a result, the mean 
magnification is not necessarily a fair indication of the magnification bias. 

Magnification bias operates differently in surveys for galaxy and cluster lenses. The 
typical separation of a galaxy lens is smaller than most optical survey resolution limits, 
and the total magnification must be used to compute the magnification bias. However, 
the total magnification is not the important quantity for cluster lensing surveys since the 
separate images are generally resolved (see also Chiba & Yoshii 1997). There are different 
approaches to making a wide separation gravitational lens survey. Firstly, candidate lenses 
may be selected from galaxy or cluster catalogues having one image of a source in good 
alignment with a potential lens. Follow up imaging can then be used to look for a second, 
fainter image. This is the approach taken by systematic surveys to date (e.g. Maoz et al. 
1997; Phillips, Browne & Wilkinson 2000). Here the magnification bias which results in 
the inclusion of the brighter image into the initial survey is important (though the overall 
bias for the survey must also consider the detectability of the faint image in follow up 
observations) . 

A second method looks for pairs of lensed images in existing data. In this instance, 
the magnification bias which also includes the fainter image (the second brightest for three 
image lenses, i.e. image 3), and thus both images into the survey is the relevant quantity. 

To compute magnification bias, a number magnitude relation for the sources is required. 
We use the function described by Kochanek (1996), based on data from Boyle, Shanks & 
Peterson (1988) and Hartwick & Schade (1990). He uses a broken-power-law form (Turner 
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1990). To simplify our calculations and to keep the conclusions as general as possible, we 
assume that the survey depth is deeper than the break B-magnitude at z s = 3 of tub ~ 19. 
The magnitude-number relation may then be written as 

dN 

— = N o l0 a ™ m (31) 
dm 

where N Q is a normalizing constant and a m is the logarithmic slope. For a survey with a 
limit m,i im , the magnification bias is defined as (Fukugita & Turner 1991) 

_ / °°f iV(< m Hm + llogM)d t i 
N{<m lim ) 

In application to a real lens survey, the magnification bias may be written more conveniently 

as 

R _ /o°° gT (< m Hm + llogM) dp, 

N(<m lim ) ' { } 

where the limits \x m %n and /z max are the range of the bright image magnification such 

that at least one other image is visible. Note that only the magnification distribution for 

the bright image is directly relevant. In this paper we adopt the simpler and non-specific 

case of magnification bias for bright and faint images based on Eqn. [32[ 

Fig. P shows the magnification distributions for the bright (left hand panel) and faint 
(right hand panel) images for Zhao profiles. The light lines show the distributions for a SIS: 

dP SIS _ 2 

dfJ,f (flf + l) 3 

and 



forfi f >0 (34) 



dP SIS 2 
d/ib (A*6 - 1) 

for comparison. The image magnification distributions for the Zhao profiles exhibit similar 
behavior to the SIS, and there is a high probability that the two images will have different 
magnifications. However unlike the SIS, the Zhao profiles have finite as 9 — > 0. The 
faint image therefore has a finite minimum magnification. 
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The magnification bias described by Eqn. 32 describes the additional number of 
lensed pairs. If the quantity of interest is the density of lensed images on the sky, then 
an additional factor of must be inserted under the integral to account for the smaller 
region of source plane sampled per unit area of lens plane. Inserting Eqn. |3l] into Eqn. [32] 
we find 

r°° dP § a 

B b = / — /i fe 2 dfi b (36) 
Jo anb 

for magnification bias due to bright images, and 

/•OO r] P 5 

B ' = l nr,"' ** (37) 

for magnification bias due to faint images. Note that the magnification bias is equal to the 
mean magnification for a m = 0.4. Since the high magnification tail of any magnification 
distribution has 

dP * 

*r <x " ■ (38) 



the integral in Eqns. and [h| is finite for a m < |, and in this range no low-magnitude 
cutoff need be assumed. 



Using Eqns. |34| and |35] the magnification biases for the faint and bright images of the 
SIS are 

B f s = f uV^- < 39 > 

and 

roc 9 5 

B ' =1 feTI?"? (40) 



Fig. |T2 shows the absolute magnification bias (left hand panel), and the magnification 
bias relative to the SIS (right hand panel) based on the bright image as a function of a m . 
Plots are shown for the values of /3 and C^fw m Tab. [I]. The light line in the left hand panel 
shows the magnification bias curve for the SIS. Fig. [TTJ shows the corresponding results for 
the faint image. Values of absolute and relative (to the SIS) magnification bias based on 
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the faint and bright image are presented in Tabs. 0, |3|, |] and | for a m = 0.2, 0.4, 0.6. For 
quasars Kockanek (1996) finds 0.27 ± 0.07. 

The absolute magnification biases for the bright images range up to a few, a few tens 
and a few hundred for a m = 0.2, 0.4 and 0.6 respectively. However, as discussed at the 
top of this section, the quantity of interest is the relative magnification bias. The values 
in Tab. ^| demonstrate the reduced amplitude of the relative (to the SIS) bias, which is 
similar for a m = 0.2, but reduced by a factors of ~ 3 and ~ 5 for a m = 0.4 and 0.6. Similar 
behavior is seen in magnification bias of faint images. The absolute magnification biases for 
the faint images produced by Zhao profiles are slightly reduced from those computed for 
the bright image. However the magnification biases for the faint images relative to the SIS 
are increased by a factor of two over the corresponding bright image case. This apparent 
contradiction arises because the faint image produced by a SIS may be arbitrarily faint. 

The magnification bias results in additional sources being available for lensing. Since 
the relative magnification bias is greater than 1 where the ratio of cross-sections is less than 
one, the magnification bias counteracts the dramatic dependence of r on (3 and Cnfw seen 
in Fig. [|. The relevant quantity is therefore the optical depth to multiple imaging r(z s )Bb 
or r(z s )Bf. Tabs. || and [7] show values of the optical depth relative to the SIS for biases 
computed using bright and faint images. Values are shown assuming a m = 0.2, 0.4 and 0.6. 
The magnification bias reduces the difference between cross-sections for different profiles, 
though the variation with (3 and Cnfw is still very dramatic. However, profiles with no 
cross-section to multiple imaging are not subject to magnification bias. Magnification bias 
therefore increases the dependence of the cross-section to multiple imaging on f3 and Cnfw 
near the multiple-imaging / no-multiple-imaging boundary. 
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7. IMPLICATIONS FOR CDM AND SIDM PROFILES 

Miralda-Escude (2000) argued that strong lensing was a powerful test of SIDM 
models. Assuming a cross-section independent of velocity, he pointed out that if SIDM 
were responsible for the cores in dwarf galaxies, then scaling arguments lead to clusters 
cores that are larger than observed. While Miralda-Escude used a scaling law that is 
not a good fit to the SIDM simulations (Dave, Spergel, Steinhardt & Wandelt 2000), his 
basic argument, that lensing is a powerful test of the model remains valid. Furthermore, 
Miralda-Escude (2000) points out that the mis-alignment of the images in the cluster lens 
MS2137-23 cannot be reconciled with a spherical lens implied by SIDM halos. 

Recently, Yoshida, Springel, White & Tormen (2000) have produced simulations of a 
cluster with M 200 = 7.4 x 1O 14 M . They model the cluster using collision-less cold dark 
matter (labeled SI), as well as weakly self interacting cold dark matter having collision 
cross-sections of a* = 0.1, 1.0 and lO.Ocm^m -1 (labeled SlW-a, SlW-b, and SlW-c). Halo 
density profiles are presented at z = for all 4 realizations. In addition, they present halo 
profiles for the simulation SlW-c at a range of different epochs. During the evolution of 
the halo, the central densities are most extreme when the cosmological expansion factor 
a (normalized to its present value) takes the values a = 0.78 and a = 1.72. The 6 halo 
profiles are reproduced in Fig. |T2| (dots). To estimate the lensing rate due to these different 
profiles we have fitted each with the Zhao profile (Eqn. |TJ) . The fitting was achieved via 
minimization of x 2 (solid lines) and Ap max (dashed lines). The constraint that (5 > 0.2 was 
imposed and the halos SlW-c (a=0.78) and SlW-c (a=1.72) were fitted assuming an r 2 oo 
(and hence M200) corresponding to that at a—1. Yoshida et al. (2000) note that the halo 
approximately doubled its M 2 oo between a = 0.73 and a = 1.72. If half this gain was made 
before (after) a = 1, then C NFW for the halo SlW-c (a=0.78) (SlW-c (a=1.72)) should be 
revised downwards (upwards) by ~ 10%. Both the fits and fractional residuals are plotted 
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in Fig. 12. The Zhao profile provides an adequate fit to halos SI and SlW-a. However for 
the larger interaction cross-sections, the core density is overestimated by the fitted profile. 
The corresponding lensing rate may therefore also be overestimated in these cases. 

The estimates for the profile parameters Cnfw an d /3 of the Yoshida, Springel, White 
& Tormen (2000) halos are plotted over contours of \ in Fig. [13|. Each parameter set 
is labeled and the upper and lower plots show points corresponding to the fits obtained 
by minimization of x 2 an d Ap max respectively. From Fig. [13] we find that r falls off very 
rapidly for small j3 and C^fw-i t ne region of interest for SIDM halos. The halos SlW-b and 
SlW-c at the current epoch do not produce any multiple images. The halo SlW-a (at the 
current epoch) lenses at a significantly reduced rate with respect to the CDM halo SI. The 
variation in the density profile at different epochs is due to successive mergers (Yoshida 
et al. 2000). The spread in parameters obtained from profile fitting is representative of 
the range of lensing strengths expected. A population of SlW-c (a = 1.72) halos would 
yield a finite lensing rate, comparable with that of SlW-a halos, while at the other epochs 
considered (a = 0.78 and a = 1.0) the halos would produce no multiple images. 



For comparison, we have also included in Fig. [13], points corresponding to parameters 
found by two complimentary studies of CDM halos. Jing & Suto (2000) have estimated of 
the spread in the distribution of concentration parameter for CDM halos by fitting the Zhao 
form to 4 cluster mass halos (M 2 oo ~ 3 — 5 x 1O 14 M ) using both (3=1 and /3 = 1.5. Their 
results are plotted in the lower panel of Fig. [TB| (labeled 1-4 as per Jing & Suto (2000)). This 
range is in agreement with the 68% range for Cnfw (extrapolated to M 2 oo ~ 1O 15 M ) found 
by Bullock et al. (2000) (plotted in the upper panel of Fig. |13|). The halo SI has a Cnfw 
which is consistent with the range of values found by these authors. It is important to note 
that the range of parameters describe halos that produce lensing cross-sections differing by 
several orders of magnitude. The lensing statistics will therefore be dominated by the more 
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concentrated halos, indicating that the medians of profile parameter distributions may not 
provide a reasonable estimation of the lensing rate for a particular dark matter model. We 
return to this point in the following section. 

More detailed modeling including evolution, non-spherical mass distributions, and a 
range of formation red-shifts, as well as observational details including magnification bias 
etc. (e.g. Turner, Ostriker & Gott 1984; Kochanek 1991) appropriate for surveys such as 
that being undertaken by the SDSS collaboration (e.g. York et al. 2000), is required before 
quantitative constraints on SIDM models can be derived. In addition, since lensing probes 
the central regions of the cluster, the contribution to the potential of a centrally dominant 
galaxy may be significant and must also be considered. We offer a preliminary discussion of 
the contribution of central galaxies in Sec. ^, however, we conclude that the sensitivity of r 
to j3 and C^fw indicates that cluster lensing statistics will provide a powerful probe of the 
self-interaction cross-section. 

While lensing constrains the dark matter properties for dark matter collisions with 
velocities of v ~ 10 3 km sec -1 , observations of dwarfs constrain dark matter properties 
with v ~ 30 km sec' 1 . If the dark matter particle is a composite particle, analogous to 
a nucleon with a low energy resonance, then its cross-section scales as a*v = const. For 
fixed cross-sections on the dwarf scale, a cross-section that is inversely proportional to 
the characteristic velocity will produce significantly smaller, less spherical cores in clusters 
and large cores in dwarfs (Yoshida, Springel, White & Tormen 2000), consistent with 
observations of cluster cores like CL0024+1654 and CL0016+16 (Avila-Reese, Firmani, 
D'Onghia & Hernandez 2000; Firmani et al. 2000). Alternatively, if the dark matter can 
self annihilate at late times (Kaplinghat, Knox & Turner 2000) then a*v = const., this 
leads to low density tri-axial cores and differing lensing predictions. 
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8. DEPENDENCE OF OPTICAL DEPTH ON THE HALO PROFILE 

DISTRIBUTION 

Recent studies of N-body simulations have described the distribution of halo 
concentrations (Bullock et al. 2000; Jing & Suto 2000; Dave, Spergel, Steinhardt & 
Wandelt 2000). Moreover, the simulations of Yoshida, Springel, White & Tormen (2000) 
demonstrate large variation in the concentration of a SIDM halo during its evolution. An 
important quantity which has hitherto been neglected in parametric cluster lens analyses is 
the distribution of halo profiles. 

The upper panels of Fig. [14| show cross-sections T ^*(l \ (thick light lines) and optical 
depths ^fe Tsistl ) (thick dark lines) to multiple imaging with respect to an SIS as a 
function of Cnfw for (3 = 1.0 (left hand panel) and (3 = 0.5 (right hand panel). We have 
assumed a m = 0.2, appropriate for the faint end of the quasar luminosity function. The 
light curves are sections through Fig. £|. The corresponding curves for ^r^^~f^j show 
the decreased slope of -Mrs T ^ n with Cnfw f° r large Cnfw, an d the increased slope 
near the no-multiple-imaging cutoff discussed at the end of Sec. [| The upper panels also 
show hypothetical Cnfw distributions having means of (Cnfw) = 6.1 (thin solid lines) 
and (Cnfw) — 10.9 (thin dashed lines). In each case one delta function and one Gaussian 
distribution (half-width of <7c NFW = 2) ar e plotted. The lower panels show the dependence 
on o~c NFW of convolutions of T ^^ z (light lines) and -^j^ T ^^ S (l ^ (dark lines) with Gaussian 
distributions of Cnfw- The solid and dashed lines correspond to convolutions with Cnfw 
distributions having (Cnfw) — 6.1 and 10.9 respectively. The convolutions in the lower left 
and right panels assumed f3 — 1.0 and (3 = 0.5. 

The optical depth shows significant dependence on the range of Cnfw- F° r example, 
the simulations of Bullock et al. (2000) yield {3=1, (Cnfw) ~ 6 and (?c NFW ~ 2 for 
CDM halos. The spread of &c NFW ~ 2 results in a 2 fold increase in optical depth to 
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multiple imaging. Much more dramatic examples are obtained where (C^fw) lies close to, 
or below the no-multiple-imaging boundary, and the significant portion of optical depth 
is produced by outliers in the concentration population. This is the region of parameter 
space important for SIDM halos. We use the SlW-b halo of Yoshida, Springel, White & 
Tormen (2000) as an example. Using the C^fw ~ ft degeneracy lines in Fig. ||, we find 
that this halo may be approximated by a Zhao profile with (3 = 0.5 and C^fw ~ 6. These 
parameters place the SlW-b halo in the no- multiple imaging regime. However, if halos 
have concentrations distributed about Cnfw ~ 6.1, the optical depth is finite. Taking the 
example of (Cnfw) — 6, we find that 2 orders of magnitude separate the optical depth due 
to distributions of halos with crc NFW = 1 an d °c NFW = 2, and that 3 orders of magnitude 
separate the optical depth between distributions having <j Cnfw = 1 and (Jc NFW = 4- Note 
that the SlW-c halo of Yoshida et al. (2000) has a concentration that differs by up to 
ACnfw ~ 8 from the present day value during its evolution. These results demonstrate 
that knowledge of the distribution of halo concentrations is an important consideration for 
cluster lensing studies. 



9. THE CONTRIBUTION OF CLUSTER GALAXIES 

The statistics of multiple imaging due to clusters of galaxies are sensitive to the inner 
most regions of the mass distribution, where the contribution of a central galaxy may be 
important. One expects that this is particularly true where the halo is less concentrated 
or or has a shallower cusp, as is the case for SIDM halos. The presence of an isothermal 
central galaxy renders the cross-section to multiple imaging finite for all halos. 

Williams, Navarro & Bartelmann (1999) have investigated the effect of including a 
massive central galaxy. They find that arc properties of the most massive clusters are 
consistent with predictions of the NFW profile, but that observed separations are too large 
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in the lower mass cluster lenses. They demonstrate that a centrally dominant galaxy with 
a dark halo can resolve the discrepancy. However there are several caveats. As noted by 
Williams et al. (1999), their cluster sample over represents the expected fraction of massive 
clusters (based on the algorithm of Press & Schechter 1974), possibly indicating a selection 
bias towards large splittings and highly elongated arcs. Furthermore, the higher resolution 
of simulations of Moore et al. (1998) show a steep core p{r) ~ r~i, indicating that the 
concentrations found by NFW may be a factor of ~ 2 too low. This reduces the required 
mass of the centrally dominant galaxy. More importantly, Bullock et al. (2000) find 
significant scatter in concentration for a halo of given mass. High concentration halos are 
more likely to lens, produce larger image splittings, and so are more likely to be observed 
as cluster lenses. In the previous section this was shown to be a significant effect and 
might give the appearance that the profiles found from N-body simulations cannot account 
for observations of cluster lenses. On the other hand, Bullock et al. (2000) also find a 
concentration that decreases faster with redshift than originally thought (Navarro, Frenk & 
White 1997). This will exasperate the problem pointed out by Williams et al. (1999) for 
cluster lenses above a redshift of ~ 1. 

On the other hand, a cluster need not have a single centrally dominant galaxy (e.g. 
CL0024+1654). Flores, Mailer & Primack (2000) have analyzed the inclusion of cluster 
galaxies on the statistics of giant arc properties. They find that observationally based 
constraints suggest that there are not enough massive galaxies in a cluster to significantly 
alter arc statistics. In a complementary study, Menghetti et al. (2000) artificially added 
galaxies to N-body simulations. When the properties of short arcs are excluded, their 
analysis suggests that the contribution of cluster galaxies is negligible. 

A proper account of the cluster density profile in the presence of a centrally dominant 
galaxy requires a full treatment of galaxy formation at the center of the cluster. The studies 
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mentioned artificially add intrinsically or extrinsically truncated parametric galactic halos 
to the cluster halo. Williams, Navarro & Bartelmann (1999) (following Navarro, Frenk 
& White (1997)) take the additional step of modifying the cluster halo adiabatically in 
response to the formed central galaxy. Profiles with large and small (3 describe CDM and 
SI DM halos respectively. The response of the dark matter to the formation of a central 
galaxy cannot be commonly approximated in these two cases. 

In CDM halos, the central galaxy may form adiabatically (see Blumenthal, Faber, 
Flores & Primack 1984), resulting in additional contraction of the cluster dark matter halo. 
In this scenario, it is possible that the central galaxy may retain its own halo if it is not 
tidally striped during formation. Note however that N-body simulations of cluster halos are 
expected to properly describe the halo of the central galaxy, since the cusp in the cluster 
density profile is formed from the infall of smaller halos. 

To reconcile observations with CDM NFW halo models, Williams et al. (1999) find 
that an additional central cluster galaxy of mass ~3x 10 12 M Q is required. With a stellar 
velocity dispersion of a* = 300 km sec' 1 , this implies that the central galaxy must retain 
its halo during formation. To find an approximate contribution to the multiple imaging 
cross-section we have taken the most simple approach and superimposed the dark matter 
halo profile of Brainerd, Blandford & Smail (1996), 

P( r ) = a n°2?l , 2V ( 41 ) 

onto the center of the Zhao profile. Here a* is the central velocity dispersion and r Q is the 
characteristic outer scale on the center of the Zhao profile. The profile has a total mass 

M 9 = (42) 

We have calculated cross-sections to multiple imaging for Zhao profiles having 

M 20 o ~ 1O 15 M in the presence of an additional central galaxy of mass M g = 3 x 10 12 M Q 
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with a central stellar velocity dispersion of = 300 km sec -1 . The resulting cross-section 
calculations are applicable to Zhao profiles with large 8 and C^fw appropriate for CDM 
halos. The left hand panels of Tabs. |8| and § show values of the cross-section to multiple 
imaging Td+h \ Zs \ and the fractional increase in cross-section Td+b ^ over a small grid of 8 

& & TSIS(Zs)' T{z a ) & h 

and Cnfw- 

In the self interacting scenario, the response of the halo to the formation of the central 
galaxy should be larger. The contraction heats up the dark matter, after which interactions 
transfer the heat outwards, resulting in cooling and additional halo contraction. An 
isothermal approximation is therefore appropriate for the response of the SIDM halo to the 
formation of a central galaxy. Consider an initial dark matter density profile resulting from 
an isothermal distribution function: 

p d (r)cxexp^ (43) 

where k is Boltzmann's constant and T the temperature. The ratio of the dark matter 
density before and after the addition of baryons is 

Pd+b §d - $d + b 2a, 

= exp — exp . (44) 

Pd Kl CT C str 

Thus the density increase of the cluster halo with a cstr ~ 1000 km sec -1 and central galaxy 
velocity dispersion cr* ~ 300 km sec -1 , the density increase is of order 20%. This density 
increase is negligible with respect to the orders of magnitude which separate the densities 
of different profiles in the inner most regions (see Fig. [I]). 

If the cluster has a SIDM halo, the central galaxy cannot have retained its own halo 
during formation. The additional galactic mass should therefore consist only of stars, and is 
thus more confined. We have calculated cross-sections to multiple imaging for Zhao profiles 
having M200 ~ 1O 15 M in the presence of an additional central of mass M g = 3 x 10 11 M Q 
(appropriate for the mass in stars of a typical central galaxy) with a central stellar velocity 
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dispersion of a* = 300 km sec -1 . The calculation is applicable to all the Zhao profiles 
assuming that the central dark-matter distribution is described solely by the Zhao profile. 
The right hand panels of Tabs. || and || show the resulting values of the cross-section to 
multiple imaging Td+b ^ Zs \ and the fractional increase in cross-section Td + b( - z ^ where a central 

& & rsisizs) 1 t{z s ) 

galaxy is included over a small grid of (3 and Cnfw- ^(z^ * s °^ or der 1 for all but the 
shallowest halo profiles. 

The addition of the extra mass increases the cross-section of all profiles with 
respect to the SIS. We have not re-normalized the mass of the Zhao-profile-plus-galaxy 
as we are interested in the relative rather than the absolute cross-sections of different 
Zhao-profile-plus-galaxies. The addition of the central isothermal galaxy sets a lower limit 
on the cross-section, and reduces the dependence of multiple imaging on the cluster profile 
parameters, though the differences are still orders of magnitude. The values in Tab. § 
may be compared directly with the magnification bias (Tabs. and f|). These show that 
the contribution of the central galaxy to the optical depth is not as important as that of 
magnification bias. Note that while magnification bias and the addition of a central galaxy 
both preferentially increase the cross-section of less concentrated clusters, the addition of 
a central galaxy reduces the magnification bias due to the increased cross-section. Thus, 
while a self-consistent calculation is yet to be done, we conclude that centrally dominant 
galaxies will not prevent lens statistics from being used as a probe of SIDM halos. 



10. CONCLUSIONS 

We have calculated the differential and total optical depths to multiple imaging, the 
average image splitting and the total magnification for a constant co-moving number 
density of unevolving generalized NFW (Zhao) profile cluster mass gravitational lenses. We 
find that the number of expected strongly lensed quasars is a very sensitive function of the 
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profile parameters. Profiles whose central density is low either due to a shallow central 
cusp, or to a scale radius that is a reasonable fraction of the virial radius have cross-sections 
to multiple imaging that are reduced by a significant factor. Moreover, the separation 
of multiple images is reduced (by a factor of a few), although the total magnification is 
significantly enhanced. We find that the resulting magnification bias does not alter our 
conclusions. 

The Zhao profile exhibits degeneracies between profile parameters with respect to 
lensing statistics. Similarly, profile fits have parameter degeneracies which are a function 
of the minimization quantity adopted. If the properties of a profile core are accurately 
reproduced by an approximately degenerate profile, then the parameter degeneracies are 
nearly equivalent to those obtained from strong lensing statistics. However, a profile that 
attempts to fit the entire halo will introduce serious uncertainty into the inferred lensing 
rates (up to an order of magnitude). This is particularly true in the region of parameter 
space where both Cnfw and (3 are small and will be an important consideration for detailed 
calculations of the lensing rate for SIDM halos based on parametric results from N-body 
simulations. 

The lensing rate is a powerful probe of SIDM in clusters of galaxies, particularly in 
the absence of a centrally dominant galaxy. We have obtained profile parameters for halos 
around clusters of galaxies by fitting Zhao profiles to the simulations of Yoshida, Springel, 
White & Tormen (2000). We find that the optical depth to multiple imaging seriously 
constrains the SIDM self-interaction cross-section. In particular, an interaction cross-section 
of 10.0cm 2 gm' 1 is rarely capable of producing multiple images during its evolution. In 
the most highly concentrated phases, a halo composed of interacting dark matter with 
this cross-section has a rate of multiple imaging 1 to 2 orders of magnitude lower than the 
corresponding typical CDM halo. Preliminary calculations show that centrally dominant 
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galaxies increase the multiple imaging cross-section, but that the increase is small with 
respect to the variation in cross-section among different profiles. The presence of centrally 
dominant galaxies should not therefore inhibit the use of the multiple imaging rate as a 
probe of SIDM. 

An important result from this study, with implications for cluster lensing studies is 
that the scatter in CDM/SIDM profile parameters obtained from N-body studies describe a 
very large range of optical depth. The lensing statistics will therefore be dominated by the 
more concentrated members of the population, rather than by the typical halo. As a result, 
an estimate of the distribution of halo profiles must be included in future studies of cluster 
lensing statistics. 
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Fig. 1. — Plots of p(r) and E(£) for Zhao profiles / (solid lines), II (dashed lines) and III 
(dotted lines). psis( r ) an d ^sis(0 f° r a singular isothermal sphere are also shown (light 
lines). The dashed lines in the lower panel show the critical density for a source at z s = 3 
and lenses at Zd = 0.5, 1.0 and 2.0. 

Fig. 2. — Plots of t(z s ) for Zhao profiles I (solid line), II (dashed line) and III (dotted 
line). Tsis( z s) for a singular isothermal sphere is also shown (light line). t(z s ) has been 
normalized by Fsis- 

Fig. 3.— Plots of * dT( * s z ' d Zd) for profiles / (solid lines), // (dashed lines) and /// (dotted 
lines) with z s = 1.5 (upper panel) and z s = 3.0 (lower panel). The light lines show the 
corresponding distributions for a SIS. 

Fig. 4. — Contour plot of r J"^] ^ with z s = 3 over a parameter space of (3 and Cnfw- The 
upper thick line is a contour of unity and marks the parameters that produce an cross-section 
equal to that of a SIS. The lower thick line marks the boundary of the region of parameter 
space where no multiple imaging is produced. The dots show the locations of profiles /, //, 
and III (left to right). 

Fig. 5. — Distributions of image separations over a grid of Cnfw =20.0 (solid lines), 12.5 
(dashed lines), 5.0 (dotted lines) and (3 =1.5 (top), 1.0 (center), 0.5 (bottom) (note that 
{Cnfw-, P) = (5.0,0.5) produces no multiple images and is not included). The light lines 
show the distribution for an isothermal sphere. 

Fig. 6. — Distributions of total image magnifications over a grid of Cnfw =20.0 (solid lines), 
12.5 (dashed lines), 5.0 (dotted lines) and (5 =1.5 (top), 1.0 (center), 0.5 (bottom) (note that 
(C NFW , (3) = (5.0,0.5) produces no multiple images and is not included). The light lines 
show the magnification distribution for an isothermal sphere. 
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Fig. 7. — A representation of the degeneracy between different Zhao profiles. Profiles with 
fixed (3 were fit to a Zhao profile (f3 = 0.5) over a range of concentration parameters 
Cnfw,o- The fitted Cnfw and average fractional difference in density for a 3-d halo are 
shown. The fractional residuals are plotted above each set of fits. The minimization 
used the following functions. Top panels: x 2 — S r 2 (p(r) — p (r)) 2 dr. Center panels: 

100 

max{\log(p(r)) - log( Po (r))|). Bottom panels: xld = S r 2 ( p(r) - p ^ r) f dr . 

100 CoV ) 

Fig. 8. — A representation of the degeneracy between different Zhao profiles. The dot- 
dashed lines are curves along which the difference between halos with different parameters 
are minimized. The minimization used the following functions. Top figure: % 2 = 
f rloo r 2 (p(r) — p (r)) 2 dr. Center figure: max(\log(p(r)) — log(p (r))|). Bottom figure: 

100 

Xrei = S r 2 ( p ^~ ( Po S r ^ ) 2 dr. The solid light lines in these figures represent the contours 

100 Po\ r ) 

of -^Kr with z s = 3. 

Fig. 9. — Distributions of bright (left hand panels) and faint (right hand panels) image 
magnifications over a grid of Cnfw =20.0 (solid lines), 12.5 (dashed lines), 5.0 (dotted lines) 
and P =1.5 (top), 1.0 (center), 0.5 (bottom) (note that (Cnfw,P) = (5.0,0.5) produces no 
multiple images and is not included). The light lines show the corresponding bright and 
faint image magnification distributions for an isothermal sphere. 

Fig. 10. — The magnification bias for bright images as a function of a m over a grid of 
Cnfw =20.0 (solid lines), 12.5 (dashed lines), 5.0 (dotted lines) and P =1.5 (top), 1.0 
(center), 0.5 (bottom) (note that (Cnfw,P) = (5.0,0.5) produces no multiple images and is 
not included). Left: The absolute bias B^. The light line shows the bias for an isothermal 
sphere B^ IS . Right: The magnification bias relative to the isothermal sphere -^h- 

Fig. 11. — The magnification bias for faint images as a function of a m over a grid of 
Cnfw =20.0 (solid lines), 12.5 (dashed lines), 5.0 (dotted lines) and P =1.5 (top), 1.0 
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(center), 0.5 (bottom) (note that (C NFW ,/3) = (5.0,0.5) produces no multiple images and is 
not included). Left: The absolute bias Bf. The light line shows the bias for an isothermal 
sphere Bf Is . Right: The magnification bias relative to the isothermal sphere 

j t> f 

Fig. 12. — Fits and residuals for the Yoshida et al. (2000) halos (dots). The solid and dashed 
lines correspond to minimization of x 2 and Ap max . Each halo is labeled as per Yoshida et 
al. (2000). 

Fig. 13. — Estimation of relative lensing rates for CDM and SIDM cluster halos. The dots 
represent the parameter combinations corresponding to fits for the Yoshida et al. (2000) 
halos (labeled, parameters in upper and lower plots used minimization of x 2 an d A/9 max ). 
The squares and triangles in the lower figure show the parameters (from Jing & Suto 2000) 
resulting from fits of Zhao profiles with (3 — 1 and (3 — 1.5 respectively to CDM halos 
(labeled 1-4 as in Jing & Suto 2000). The range bar in the upper figure shows the range 
(68%) of concentrations found for the (3 = 1 profile by Bullock et al. (2000) (extrapolated 
to M 20 o ~ 10 15 M Q ). These are superimposed on contours of T ^^l j with z s = 3. 

Fig. 14. — Upper panels: - 1 — (thick light lines) and ■^^— L - (thick dark lines) as a function 
of C NFW for (3 = 1.0 (left hand panel) and (3 = 0.5 (right hand panel). The thin solid and thin 
dashed curves show example Cnfw distributions having (Cnfw) = 6.1 and 10.9. In each 
case one delta function distribution and one Gaussian distribution (half-width of o~c NFW = 3) 
are plotted. Lower panels: Convolutions of the cross-section ta J (light lines) and optical 
depth ^Irg-^j^ (dark lines) with example Gaussian Cnfw distributions as a function of 
<Tc NFW - The solid and dashed lines correspond to convolutions with Cnfw distributions 
having (Cnfw) — 6.1 and 10.9 respectively. The Convolutions in the lower left and right 
panels assumed (3 = 1.0 and (3 = 0.5 respectively. We have assumed a m = 0.2. 
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Table 1: Values for T f z \ , (A0(z s )) and {n(z s )) over a small grid of (3 and CjVFiy- 



5.0 12.5 20.0 5.0 12.5 20.0 5.0 12.5 20.0 

0.5 0.000 0.008 0.174 NA 26.52 46.64 NA 75.41 18.46 

P 1.0 0.001 0.199 0.918 11.56 43.06 57.13 100.1 16.11 8.105 

1.5 0.226 1.715 4.263 34.91 57.79 68.71 12.33 5.732 3.861 
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Table 2: Values for the bright image magnification bias Bf, over a small grid of (3 and Cnfw- 



Cnfw Cnfw Cnfw 

5.0 12.5 20.0 5.0 12.5 20.0 5.0 12.5 20.0 



0.5 NA 5.53 2.50 NA 35.57 7.35 NA 282.1 27.45 
P 1.0 6.27 2.53 1.89 45.32 7.50 4.13 398.1 27.91 11.18 
1.5 2.32 1.68 1.45 6.24 3.22 2.38 20.50 7.54 4.62 



Table 3: Values for the bright image magnification bias relative to the isothermal sphere -Jsrz over 
a small grid of [3 and Cnfw- 



Cnfw Cnfw Cnfw 

5.0 12.5 20.0 5.0 12.5 20.0 5.0 12.5 20.0 



0.5 NA 3.76 1.70 NA 13.34 2.76 NA 49.51 4.82 
P 1.0 4.27 1.73 1.29 16.99 2.81 1.55 69.88 4.90 1.96 
1.5 1.57 1.15 0.99 2.34 1.20 0.89 3.60 1.32 0.81 
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Table 4: Values for the faint image magnification bias Bf over a small grid of and Cnfw- 



Cnfw Cnfw Cnfw 

5.0 12.5 20.0 5.0 12.5 20.0 5.0 12.5 20.0 



0.5 NA 4.94 2.35 NA 29.23 6.55 NA 227.5 24.00 
P 1.0 5.94 2.36 1.56 41.98 6.59 2.97 364.7 24.12 7.90 
1.5 1.96 2.22 0.91 4.61 1.87 1.10 14.76 4.25 2.13 



Table 5: Values for the faint image magnification bias relative to the isothermal sphere -Mrs over 

B f 

a small grid of (3 and Cnfw- 



Cnfw Cnfw Cnfw 

5.0 12.5 20.0 5.0 12.5 20.0 5.0 12.5 20.0 



0.5 NA 6.73 3.20 NA 30.90 6.93 NA 101.7 10.74 
P 1.0 8.10 3.21 2.13 43.84 6.96 3.14 163.1 10.79 3.53 
1.5 2.67 1.66 1.24 4.88 1.98 1.16 6.60 1.90 0.95 
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Table 6: Values for the bright image magnification bias weighted cross-section relative to the 
isothermal sphere j B over a small grid of (5 and Catfvk- 



Cnfw Cnfw Cnfw 

5.0 12.5 20.0 5.0 12.5 20.0 5.0 12.5 20.0 



0.5 0.000 0.030 0.296 0.000 0.107 0.480 0.000 0.396 0.839 
P 1.0 0.004 0.344 1.184 0.017 0.559 1.423 0.070 0.975 1.799 
1.5 0.355 1.972 4.220 0.529 2.058 3.794 0.814 2.264 3.453 



Table 7: Values for the faint image magnification bias weighted cross-section relative to the 
isothermal sphere gt/is over a small grid of f3 and Cnfw- 



Cnfw Cnfw Cnfw 

5.0 12.5 20.0 5.0 12.5 20.0 5.0 12.5 20.0 



0.5 0.000 0.054 0.557 0.000 0.247 1.206 0.000 0.814 1.869 
P 1.0 0.008 0.639 1.955 0.044 1.385 2.883 0.163 2.147 3.240 
1.5 0.603 2.847 5.286 1.103 3.396 4.945 1.492 3.259 4.050 



-58- 



Table 8: Values for Td + b { z °) over a small grid of 3 and Cnfw- Each cluster contains a central 

T SIS\ z s) ' ° ' vv 

galaxy of mass 3 x 1O 12 M (left) or 3 x 1O U M (right) having a central velocity dispersion of 
&CD = 300 km seer 1 . 



M g = 3x 10 12 M o M g = 3 x lO n M 



Cnfw Cnfw 
5.0 12.5 20.0 5.0 12.5 20.0 



0.5 0.008 0.026 0.233 0.008 0.015 0.188 
(3 1.0 0.014 0.266 1.054 0.010 0.217 0.952 
1.5 0.306 1.916 4.581 0.257 1.778 4.364 



Td+b{Zs) 
T(Zs) 

galaxy of mass 3 x 1O 12 M (left) or 3 x 1O 11 M (right) having a central velocity dispersion of 



Table 9: Values for ^777-) , over a small grid of 3 and Cnfw- Each cluster contains a central 



&CD = 300 km sec 1 . 



M g = 3x 1O 12 M Mg = 3 x lO n M 



Cnfw Cnfw 
5.0 12.5 20.0 5.0 12.5 20.0 



0.5 NA 3.250 1.339 NA 1.875 1.080 
P 1.0 14.00 1.337 1.148 10.00 1.090 1.037 
1.5 1.354 1.117 1.075 1.137 1.037 1.024 



